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We present for the first time solutions in the gauged U(l) x U(l) model of Witten describing vortons - 
spinning flux loops stabilized against contraction by the centrifugal force. Vortons were heuristically described 
many years ago, however, the corresponding field theory solutions were not obtained and so the stability issue 
remained open. We construct explicitly a family of stationary vortons characterized by their charge and angular 
momentum. Most of them are unstable and break in pieces when perturbed. However, thick vortons with small 
radius preserve their form in the 3 + 1 non-linear dynamical evolution. This gives the first ever evidence of 
stable vortons and impacts several branches of physics where they could potentially exist. These range from 
cosmology, where vortons could perhaps account for the cold dark matter, to condensed matter physics. 

PACS numbers: ll.lO.Lm, 11.27.+d, 98.80.Cq 



More than 25 years ago Witten introduced the idea of su- 
perconducting cosmic strings in the context of a field theory 
model that can be viewed as a sector of a Grand Unification 
Theory (GUT) [1]. The model admits classical solutions de- 
scribing strings (vortices) whose longitudinal current can at- 
tain astronomical values, implying interesting effects (see [2] 
for a review). 

Soon after, it was realized that superconducting strings 
could form loops whose current would produce an angular 
momentum supporting them against contraction [3]. If sta- 
ble, such cosmic vortons should be of considerable physical 
interest. However, until very recently it was not clear if vor- 
tons are stable or not, since the underlying field theory so- 
lutions were not known. Various approximations were used 
to describe vortons, for example, by viewing them as thin 
and large elastic rings [4]. At the same time, it was realized 
that objects similar to vortons could potentially exist also in 
other branches of physics, as for example in condensed matter 
physics [5], or in QCD [6]. Since superconducting strings ex- 
ist in the Weinberg- Salam theory [7], vortons are potentially 
possible also there. 

The first field theory solutions describing stationary vortons 
were found in the global limit of Witten 's model, when the 
gauge fields vanish [8]. These vortons have approximately 
equal radius and thickness and remind a Horn torus. Solu- 
tions describing thin and large vortons were later found as 
well, however, when perturbed, thin vortons turn out to be 
dynamically unstable and break in pieces [9] . Although dis- 
couraging, this result is actually quite natural, since thin vor- 
tons can be locally approximated by straight strings, while the 
latter are known to become unstable for large currents [2]. 

However, a more close inspection reveals that unstable 
modes of superconducting strings have a non-zero minimal 
wavelength [10], as in the case of the Plateau-Ray leigh insta- 
bility of a water jet [11]. Therefore, imposing periodic bound- 
ary conditions with a short enough period should remove all 



instabilities. As a result, thick vortons made of short string 
pieces have chances to be stable. 

In this letter we present for the first time stationary vorton 
solutions in the gauged Witten's model, and our vortons are 
thick. To study their stability, we simulate their full 3 + 1 non- 
linear dynamics in the limit of vanishing gauge couplings. We 
find that most of them are unstable, however, thick vortons 
with a large charge and the smallest possible radius are stable. 
By continuity, it follows that vortons with non-zero but small 
gauge couplings should be stable as well. 

We therefore present the first evidence of stable vortons, 
whose features turn out to be quite different from those pre- 
dicted by the effective theories. This can impact several 
branches of physics where vortons could potentially exist, 
ranging from cosmology and high energy physics to con- 
densed matter physics. 

The model of Witten.- This is a theory of two Abelian 
vectors A^"^ interacting with two complex scalars {cl = 
1,2) with the Lagrangian 

c = -\y, fI^'Jf^-^ + Y.^D^^arD^^a -V. (1) 

a a 

Here the gauge field strengths are Fjn) — O^A^ ^ — Oi^A^ \ 

the gauge covariant derivatives D^^a = [d^ + igaA^i^^)^a 
with gauge couplings Qa, and the potential is 

^ = E T - + 7l*in*2p - ^ , (2) 

a 

where r^i = 1. If 47^ > A1A2 and 27 > A2?72 then the global 

minimum of the potential (vacuum) is achieved for | = 1 

and <l>2 = 0. Fields A^^l\ <l>i, ^2 are massive with masses, 

respectively, m\ = 2g\,m\ = Xi^m^ =7—^X2^2 whereas 
(2) 

Ajji is massless and can be identified with electromagnetic 
field. The theory has a local U(l) x U(l) invariance and two 



Figure 1. (Color online) - Profiles of the stationary vorton solution for Q = 1500 and n = 1, m = 1 for the parameter values Ai = 41.1, 
(A2, 772) = (30, 1), 7 = 20, and gi — g2 — 0.01. The configuration is very compact and close to a Horn torus. The first panel displays 
constant energy surfaces. The second panel shows surfaces of constant |$i|^ and the magnetic field B*^^^ (cones). One has — 0. The 
third panel shows |$2|^ isosurfaces and the electromagnetic current ^'2 (cones). The last panel shows the electric field E^'^^ (cones), while the 
isosurfaces show its magnitude. The red color corresponds to large values and the blue color to small values. 



Noether currents, = 2Re(z^*I^^^a),with two conserved 
charges J j^^x. The Euler-Lagrange equations are 

3V 

(9^F(«)'"^=5,i:, D^D^$, + ^^$, = o. (3) 

Assuming cylindrical coordinates = (t, p, (p), we make 
the ansatz for stationary, axially symmetric fields, 

<l>i = Xi + iYi , ^2 = {X2 + iY2) exp{i{ujt + m(p)} , (4) 

where Xa, Ya as well as A^^-' depend on p, z, and we impose 
the gauge condition Ap^^ =0. Here m is an integer winding 
number and is a frequency. The fields should be globally 
regular and the energy should be finite, which requires that 
at infinity Xi ^ 1 while all other amplitudes approach zero. 
At the symmetry axis, p = 0, the amplitudes X2,l2,^^^^ 
vanish, while for the other amplitudes the normal derivative 
8/ dp vanishes. Under the reflection z — z the amplitudes 
Ya are odd whereas all the others are even. 

The choice of the ansatz implies that the first Noether 
charge vanishes, while the second one is 

Q = 2 d^x (X| + Yi) (oj - ^2Af ^) . (5) 

The energy is ^ = / T^d^x and the angular momentum 

J = j Tl&^x = mQ, (6) 

where the energy-momentum tensor is obtained by varying 
the metric tensor, = 2^^^a£/%^^ - (5^ C. In the above 
formulas all fields and coordinates are dimensionless. If 77 is 
the energy scale, then the dimensionful (boldfaced) quantities 
are = ''7^0, A^'^^ = rjA^^\ x^ = :si^r], E = rjE, hence 
T] is the asymptotic value of $1. 

Stationary vortons.- Inserting the ansatz (4) to the field 
equations (3) gives, after separating t and (p variables, an el- 
liptic system of 10 non-linear PDE's for the 10 functions of 



p, z. We solve these equations with two different numeri- 
cal methods: using the elliptic PDE solver FIDISOL based 
on the Newton-Raphson procedure [12], and also minimizing 
the energy within a finite element approach provided by the 
Freefem++ library [13]. 

The resulting vorton solutions have a toroidal symmetry, 
with non-trivial windings along both torus generators. Apart 
from the azimuthal winding number m, there is a second in- 
teger, n, determining the winding of phase of around the 
boundary of the (p^z) half-plane. If n 7^ then vanishes 
at a point (po^O) corresponding to the center of the closed 
vortex forming the vorton, and the phase of <l>i winds around 
this point. We have constructed vortons for n = 0, 1, 2 and 
m = 0, . . . , 12. Solutions for n = m = are spherically 
symmetric, they have J = and are similar to Q-balls [14]. 

The vorton can be visualized as a toroidal tube confining a 
magnetic flux of B^^"^ = V x A^^^ since ^1 « inside the 
tube and thus the first U(l) is restored. ^^^^ is non-zero inside 
the tube, giving rise to a charged condensate and to a persis- 
tent electric current along the tube. The current creates a mo- 
mentum along the azimuthal direction, which gives rise to an 
angular momentum along z-direction. Outside the vorton tube 

the massive fields A^f}\ ^1,^2 rapidly approach their vacuum 

(2) 

values and there remains only the long-range massless Ajj, ^ 
generated by the electric current confined inside the vorton 
tube. At large r = \J p^ ^ z^ one has J^^^ = Q/(47rr) + . . . 
and A^^^ = /i sin 6>/r^ + ... , therefore, from far away the vor- 
ton looks like a superposition of an electric charge Q with a 
magnetic dipole p. 

Fig. 1 shows the 3D solution profiles for an m = 1 vorton. 
One can see that the vorton tube is very thick and compact. 
The vortex magnetic field B^^"^ = V x A*^^^ and the electric 
current are tangent to the azimuthal lines. The electric field 
E^'^'> = — VAp^ is mostly oriented along the radial direction 
and supports a non-zero flux at infinity, Q = § dE^'^^ ■ dS = 
g2Q. The massless magnetic field B^'^'^ = V x A*^^^ at large r 
is of magnetic dipole type. 
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Vortons can be labeled by their charge Q and the integer 
'spin' m = J/Q (assuming that n = 1). The vorton radius 
po is not very sensitive to the value of Q but increases rapidly 
with m, so that for large m vortons are thin and large, with 
the radius much larger than the thickness. On the other hand, 
increasing Q increases the thickness of the vorton tube, so that 
for large Q vortons are thick and look almost spherical. 




Figure 2. (Color online) - Left panel shows E against uj for the val- 
ues of gauge couplings (91,92) = (0.08,0.02) (blue), (0.01,0.01) 
(red), (0.08, 0.02) (green). Right panel shows E, Q, /i against 
91 = 92 for a fixed uj = 0.804. In both cases Ai = 41.1, A2 = 40, 
7 22.3, 772 = 1, n = 1, m = 2. 

The frequency u can be used instead of Q to characterize 
the solutions, which exist only within a finite frequency range, 
uj- < uj < ujj^. Both E and Q diverge for uj ^ uj± and 
have a minimum in between, as shown in Fig. 2. Vortons can 
be viewed as boson condensates, which is why their charge 
cannot be too small, since the boson condensation is not ener- 
getically favored for small quantities of the field quanta. For 
solutions with large Q and uj close to uj- the value of E/Q 
is less than for free particles, therefore they cannot disinte- 
grate into particles. In addition, for these solutions one has 
E{m) < mE{m = 1), so that m > 1 vortons cannot dis- 
integrate into m = 1 vortons. However, the possibility of 
disintegrating into Q-balls with n = m = is not excluded. 

In the limit of vanishing gauge coupling constants ga the 
gauge fields vanish and the solutions reduce to global vortons 
obtained in [8]. These can be visualized as closed tubes in- 
side which <l>i ^ but $2 is non-zero and creates the global 
azimuthal current, while outside ^ 1 and ^2 ^0. For 
Qa 7^ the gauge fields are excited and increase the total en- 
ergy and charge, as shown in Fig. 2. 

We find that solutions do not exist for arbitrary values of the 
model parameters A^, 7^2, 7, 9a but only for some regions in 
the parameter space. For example, fixing all parameters and 
also UJ and varying gi = g2, the solutions exist only within a 
finite range of gauge couplings, as is seen in Fig. 2. 

Dynamical vortons.- The key question is the vorton sta- 
bility, and we therefore analyze their non-linear 3 + 1 dynam- 



ics. In doing this, we consider only the global vortons without 
gauge fields, since simulating dynamics of the fully gauged 
vortons would require too much computer power. However, 
we expect the results obtained in the global case to apply to the 
fully gauged vortons as well, at least for small enough gauge 
couplings ga . Indeed, if = then the vorton is made of two 
scalars with global currents j^{^a)- For 9a 7^ ^ the cur- 
rents give rise to the source gaJa in the gauge field equations 
(3), hence for small ga one will have A^^^ ~ ga. The backre- 
action of the gauge fields on the scalars will be then propor- 
tional to g'^ and can be neglected as compared to the reaction 
of the scalars on themselves. As a result, the gauged vortons 
for small ga are essentially global vortons and the gauge field 
dynamics can be neglected. This is confirmed by the numer- 
ics, as can be seen in the right panel of Fig. 2. 

The dynamical evolution of global vortons was previously 
studied in [9] for large 'spins' m, when vortons are thin and 
large. It was found that when perturbed, the vortons rapidly 
develop pinching deformations which break them in pieces. 
This is quite natural, since thin and long vortons can be lo- 
cally approximated by straight vortex pieces, but the latter are 
known to be unstable, unless they are so short that there is no 
room for the instability to settle in [10]. Therefore, to get rid 
of instabilities one has to 'make vortons' of as short vortex 
pieces as possible, in which case the vorton radius would be 
the smallest. In other words, the 'spin' m should be small. 

To verify this conjecture, we considered a hyperbolic evolu- 
tion scheme based on an implicit /3-Newmark finite difference 
approximation. The space dependence is implemented within 
a finite element method providing a triangulation of a spher- 
ical volume in 3-space. At the boundary of the volume the 
fields are held at their vacuum values, which insures the ab- 
sence of outgoing radiation (details of our numerical method 
are given in the Appendix). The initial configuration is a sta- 
tionary, axially symmetric vorton solution. This configuration 
becomes automatically perturbed by the space discretization, 
and this triggers a non-trivial temporal evolution. 

The natural timescale is set by the value of uj of the un- 
derlying vorton solution, which is of order one, and so one 
can choose the timescale to be one. We then integrate with 
the timestep At = 0.1 and find that the m > 3 vortons very 
quickly become strongly deformed and then break in pieces. 
The time they take to break decreases rapidly as the 'spin' 
grows. For example, the m = 6 solution shown in Fig. 3 
breaks already for t ~ 20, while the one with m = 3 survives 
till t ^ 4:0 (see [15] for the online videos). 

The products of the vorton decay are typically two or more 
out-spiraling fragments of spherical topology (Q-balls) carry- 
ing fractions of the initial Noether charge. The angular mo- 
mentum of the initial vorton transforms into the orbital an- 
gular momenta of the fragments. The number of fragments 
depends on the initial charge and angular momentum. 

We therefore conclude that thin and large vortons are un- 
stable, thus confirming the result of [9] . One should say that a 
different conclusion was previously made in Ref. [16], where 
thin vortons were found to be stable. Since neither our anal- 
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Figure 3. (Color online) - Snapshots of the vorton time evolu- 
tion. The first line shows a constant |$2|^ surface for the m = 6 
solution for t — 0,17,21,28. The second and third lines show 
the constant |$i |^ surfaces and the electric current, respectively, for 
m = 3 vorton for t = 0, 31, 28, 47 and for the m = 1 vorton for 
t — 0, 104, 153, 239. The third line shows the stable solution for 
Q — 6000 and for the same values of Aa,?72,7 as in Fig. 2. 



ysis nor that of [9] confirm this, it is possible that the con- 
clusion of [16] is an artifact of modifying the scalar potential 
made in that work in order to improve the stability behavior. 
The authors of [16] also admit that their results are extremely 
sensitive to the initial conditions (unlike in our case), hence it 
is possible that they find some long-living metastable states. 

We finally turn to vortons with m = 1 , 2 and choose a large 
value Q, in which case the vortons are compact and thick. 
For m = 2 we cannot make a definite conclusion, since the 
vortons do not actually break but sometimes become strongly 
deformed. Although each time they manage to restore their 
original shape, it is not excluded that they will break in the 
long run. 

However, nothing at all happens to the m = 1 vortons. As 
the time goes, they only move slowly in the box, sometimes 
reflecting from the boundary, without any visible change of 
shape. We integrated up to t ~ 10^ (which requires weeks 
of runtime) without noticing any change in their behavior. 
We also checked that increasing the size of the box does not 
change anything, so that one cannot say that the boundary has 
a stabilizing effect. We therefore conclude that vortons with 
the lowest 'spin' and a large enough charge are stable. Intu- 
itively, this is because they are so thick that they are 'hard to 
pinch'. 

Interestingly, a very similar conclusion was made for the 
'spinning light bullets' which share many properties with vor- 
tons [17]. These are non-relativistic solutions for a complex 
scalar field with a (/:?-dependent phase (like for $2), and they 
also have toroidal profiles. The profile of their energy E{uj) 
is similar to that shown in Fig. 2. It was found that the m = 1 
solutions with a large charge and uj close to uj- preserve their 



shape in the 3 + 1 dynamical evolution [17]. Exactly the same 
statement applies for our relativistic vortons. 

In summary, for the first time since the vortons were heuris- 
tically described almost 25 year ago [3], we present the under- 
lying stable solutions within the U(l) xU(l) gauge field theory 
of Witten [1]. These vortons are quite different from what was 
expected before - instead of being thin and large rings, they 
look like thick doughnuts of almost spherical shape. 

We can now make some estimates. Assuming the original 
motivation of Witten, the energy scale should be of the GUT 
magnitude, r] ^ 10^"^ GeV. Using the value E ^ b x 10^ for 
estimates (see Fig. 2), it follows that vortons are very small 
and extremely heavy, with E~5xl0^^ GeV, which is not far 
from the Planck energy. If some versions of GUTs containing 
a U(l) X U(l) sector indeed applied in the past and the vortons 
have survived till present, then they could perhaps account for 
the non-relativistic cold dark matter. Since ^2 can be small or 
vanish, the vorton electric field is not necessarily very large. 

For = and stationary fields Eq.(3) can formally be 
interpreted as the non-relativistic Gross-Pitaevskii equation 
for a two-component Bose-Einstein condensate (BEG) [8]. 
Therefore, our solutions describe also vortex loops in BECs. 
Equally, Eq.(3) can be viewed as describing a condensate 
of K-mesons in QCD, hence our solutions describe the K- 
vortons, whose existence was conjectured in [6]. Setting the 
scale to be ~ 200 MeV gives for their energy E ~ 1 TeV. 
Although in each particular case one has to adjust the parame- 
ters to the physical values, the conclusion is that the existence 
of stable vortons impacts several branches of physics. 
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APPENDIX. NUMERICAL DETAILS 

We used three completely different numerical techniques 
for the vorton construction. First, we performed our calcu- 
lations using the elliptic PDE solver FIDISOL based on the 
iterative Newton-Raphson method [12] within a finite differ- 
ence scheme. This method was used to solve the equations (3) 
for the stationary, axially symmetric fields (4). Next, we re- 



produced the same solutions using the energy minimization 
scheme within a finite element framework provided by the 
Freefem++ library [13]. Finally, we simulated the hyperbolic 
evolution of the full 3 + 1 theory, also using the Freefem++ 
library [13]. Below we shall describe the two last methods. 

Energy minimization 

The energy-momentum tensor of Witten's model (1) is 

a a 

+ Y,{DM*D'''^a-StiC, (A.l) 

a 

and the energy E = J T/d^x. Using the axially symmetric 
ansatz in the radial gauge, 

A^^^dx^ = V^'^^dt + W^^'^dip + Z^^'^dz , (A.2) 
^1 = Xi + iYi , ^2 = {X2 + iY2) exp{i(cjt + mcp)} , 

where X^, Fa, V^^'K W^^'K Z^^^ depend on p, z, gives 



+ ((c. - g2V^'^y + ^ - 92W^"^y + {d2Z^"^y + Y {^2 + Yi - 24)^ (X| + r,^)} d^x , (A.3) 



where V = (dp^dz) denotes the gradient operator in the 
(p, ^) -plane. We find stationary solutions by minimizing the 
energy for a fixed Noether charge. To this end, we express the 
frequency uj in terms of Q via rewriting Eq. (5) in the form 



with 



so that 



^1 = J (- 
^2 = I d'xg2V^^^ (Xl^Yi) 



Q + 2E2 
2Si 



(A.4) 



(A.5) 



(A.6) 



which should be inserted to (A.3). The energy is then mini- 
mized at a fixed Q. 



Since the integrand in (A.3) does not depend on (p and is 
invariant under z ^ —z,we choose the integration domain to 
be a quadrant of radius r^^^, 

^ := e [0,rj X [0,rj Vp^^^ < r^,,| , 

(A.7) 

this is illustrated in Fig. 4. Regularity conditions at the sym- 
metry axis require that at p = one has 

X2 = 1^2 = 0, dpXi=dpYi=0, 

W^""^ = , V^^) = SpZ^^) = . (A.8) 

The reflexion symmetry in the equatorial plane implies that at 

z = 



dzXa=0, Ya=0, 



(A.9) 



Since the massive fields approach their asymptotic values ex- 
ponentially fast, we assume that for ^J ^ z^ — r^^^ the 
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fields take on the vacuum values 



dynamical equations for the scalars read in the global limit 



Xi = 1 , Yi = X2 = I2 = , and A^^^ = 0, (A.IO) 

but then we refine this by imposing mixed boundary condi- 

(2) 

tions for the massless Af^ ^ to take into account its slow falloff . 




Figure 4. 
domain. 



(Color online) - Triangulation over the computational 



The variational problem of minimizing (A. 3) is defined for 
numerical computations by adapting a finite element formula- 
tion provided by the Freefem++ library [13]. The discretiza- 
tion within this formulation is done via a (homogeneous) tri- 
angulation over ft, based on the Delaunay-Voronoi algorithm 
(see Fig. 4). The functions are expanded with respect to a 
continuous piecewise quadratic polynomial basis on each tri- 
angle. The accuracy of the method is controlled by the num- 
ber of triangles, by the order of expansion in the basis for each 
triangle, and also by the order of the quadrature formula for 
the integral over the triangles. Once the problem is mathemat- 
ically posed, a numerical optimization algorithm is used to 
solve the nonlinear variational problem of finding the minima 
of E. We used a nonlinear conjugate gradient method. The 
algorithm is iterated until the relative variation of the norm of 
the gradient of the functional E with respect to all degrees of 
freedom is less than 10~^. Virial relations which should be 
fulfilled by the solutions are then checked. 



Hyperbolic evolution 

To investigate the dynamical stability of the global vortons, 
the axially symmetric solutions are used as the initial data for 
a 3 + 1 evolution code. This time the integration domain Q 
is an (open) bounded subset of with the boundary and 
with the attached normal derivative dn- We choose Q to be 
the interior of a sphere of radius r^^^. The discrete version is 
obtained by a homogeneous triangulation over ft. The follow- 
ing numerical scheme was used to evolve the vortons. The 



dV 



(A. 11) 



Here A = is the Laplace operator and V the gradient 
operator in three space dimensions in Cartesian coordinates. 
The weak formulation of the hyperbolic problem (A. 1 1) is ob- 
tained by multiplying the equation by test functions Wa, inte- 
grating over Q and applying the Stokes formula, 

Jn Jn Jn d\^a\^ 

Wadn^a=0. (A.12) 



/ 

JdQ 



The time discretization is achieved using a /3-Newmark 
scheme, which converts (A.12) to 



/ 



Wa 



^[n +l]_^^[n]^^[n-l] 

dt2 



/ Wadn^^:^=^, (A.13) 



where V^'^'^ are the value of ^a^V at the time moment 
to + ndt. Here < /3 < 1 and the scheme is unconditionally 
stable for /3 > 1/4. Typically we choose /3 = 1/4, which cor- 
responds to the constant average acceleration method, while 
choosing /3 = 1/2 would reproduce the Crank-Nicolson 
scheme. We assume that dt^a = ^ at dVt, so that the bound- 
ary values of fields are frozen. 

The time-discretized equations (A.13) can be rewritten as a 
recurrence 

^[n+l] ^ ^-1 ^^^[n] + C) - (A.14) 

where v\'^'^ = {i^-"^^} is the vector containing all discretized 
degrees of freedom at the n-th time step. A, B are matrix 
bilinear operators and C is a vector containing the nonlinear- 
ities. 



I (2^ 

JQ 



on 



(A. 15) 



with Pi = p&t^ and /^s = (1 - 2p)&t^. Using the station- 
ary, axially symmetric solutions obtained by minimizing the 
energy, the recurrence is initialized by 



^\=XMz\)+isg{z)Yi{p,\z\), 
$0 = (X2(p, \z\) + isg{z)Y2{p, \z\)y"'^ , $^ = ^It 



$1 = $0 

1 ^O^icjdt 

(A. 16) 
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The time step dt is typically chosen to be 0.1 and the con- figuration is evolved for several hundreds of internal periods 

T = 27r/uj. 



